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ABSTRACT 

We present the first determination of the radial velocities and metallicities of 78 red giant stars in the 
isolated dwarf irregular galaxy WLM. Observations of the calcium II triplet in these stars were made 
with F0RS2 at t he VLT-UT2 in two separa ted fields of view in WLM, and the [Fe/H] values were 
conformed to the ICarretta fc GrattonI (|1997| ) ([Fe/Hjccgy) metallicity scale. The mean metallicity is 
([Fe/H]) = — 1.27±0.04 dex, with a standard deviation of cr = 0.37. We find that the stars in the inner 
field are more metal rich by A[Fe/H]= 0.30 ± 0.06 dex. These results are in agreement with previous 
photometric studies that found a radial population gradient, as well as the expectation of higher 
metallicities in the central star forming regions. Age estimates using Victoria-Regina stellar models 
show that the youngest stars in the sample (< 6 Gyr) are more metal rich by A[Fe/H]= 0.32 ± 0.08 
dex. These stars also show a lower velocity dispersion at all elliptical radii compared to the metal-poor 
stars. Kinematics for the whole red giant sample suggest a velocity gradient approximately half that 
of the gas rotation curve, with the stellar component occupying a thicker disk decoupled from the 
HI rotation plane. Taken together, the kinematics, metallicities, and ages in our sample suggest a 
young metal-rich, and kinematically cold stellar population in the central gas-rich regions of WLM, 
surrounded by a separate dynamically hot halo of older, metal poor stars. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: dwarf — galaxies: individ- 
ual (WLM) — galaxies: kinematics and dynamics 
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1. INTRODUCTION 

Dwarf galaxies play a critical role in our understand- 
ing of the assembly of galaxies in LCDM cosmolo- 
gies. With masses of 10^ to 10^ Mq, these galax- 
ies are thought to be similar to the proto-galactic 
fragments tha t merged and colla. p sed to forni large 
galaxies (e.g . iNavarro et all Il997t iMoore et all Il999l : 
iMadau et al.l l2001h . Analysing the survival of these 
low mass obje cts, particularly through reionisation 
()Ricotti fc Gned in 2005; Gnedin & Kravtsov 2006), is 
crucial to constraining galaxy formation models. For 
example, what was the minimum halo mass that could 
retain its baryons through reionisation? Theoretical con- 
straints on low mass galaxies are also provided by exam- 
ining the detailed characteristics of the various types of 
dwarf galaxies (irregulars, spheroidals, and the new ul- 
tra faint dwarfs) and the connections between them. Are 
dwarf irregular galaxies simply dSph which have under- 
gone recent gas mergers (|Demers et al.ll2006l:rBrook et al.l 
|2007( )? Are transition galaxies gas-rich dwarf galaxies 
being subject to ram pressure strippin g (e.g., Pegasus 
dwarf galaxv: iMcConnachie et al.ll2007f )? Are signatures 
of thick stellar disks or spheroidal components expected 
in dwarf irregular galaxies? How will the dynamics of the 
stellar populations compare to the gas motions in the low 
mass galaxies? This line of research can be carried out on 
the nearby Local Group galaxies, thus defining near-field 
cosmology. 

Dwarf irregular (dirr) galaxies hold a special status in 



the analysis of the Local Group galaxies because most are 
relatively isolated. Detailed studies of the nearby dwarf 
spheroidal galaxies have revealed complex and varied star 
formation his tories that have left behind distinct stellar 
populations (iTolstov et al.l 12004 iBattaglia et"aI1 120061 : 
iBosler et al.ll2007l) . However, interpretation of the kine- 
matics of these stellar populations, and therefore evolu- 
tion of the nearby dSph galaxies, is complicated by the 
fact that they exist in the dark matter halos of the MW 
and M31 (i.e., their stellar populations have likely been 
tidally stirred) . On the contrary, dwarf irregular galaxies 
are relatively isolated low mass galaxies. Evolved stellar 
populations in dirrs may prove to be excellent tracers of 
the dynamical history of low mass dwarf galaxies at early 
times, and therefore excellent comparisons for galaxy for- 
mation models. 

However, there have been few studies of the kinemat- 
ics of the stellar populations in isolated dIrrs and no de- 
tailed spectroscopic analysis of their older stellar pop- 
ulations due to their distances.^ This is significant be- 
cause various galaxy formation scenarios predict different 
characteristics for the stellar populations in e arly dwarf 
galaxies; e.g., simulations bv lMaver et al.l (2006) predict 
disk like systems that become more spheroidal through 
tidal interaction s and ram pressure stripping, whereas 
iKaufmann et al.l ()2007l ) suggest that dwarf galaxies start 
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^ We note that ITolstov et al.l l|2001l ') examined the calcium II 
triplet feature in 23 RGB stars in the closest dIrr, NGC 6822 
(^TJJGS ~ 21), and found most stars were young and metal-rich, 
representative of the dominant young population at the brightest 
magnitudes. 
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out as thick, puffy systems and through gas losses and 
tidal interactions become more disk like. 

In this paper, we present the first spectroscopic analy- 
sis of the calcium II triplet (CaT) feature in a sample of 
RGB stars in the dirr galaxy WLM. WLM is a typical 
low-luminosity, high gas fraction late-type dwarf irregu- 
lar galaxy. A summary of its fundamental parameters 
is listed in Table 1. The nearest neighb our to WLM is 
the C etus dSph, which lies 200 kpc away (jWhiting et al.l 
I1999D . WLM's dista nce to the Milky W ay is estimated 
at - 932 ± 33 kpc (jMcConnachie et al.l |2005,), and its 
separation from M31 is ~820 kpc, therefore this is one 
of the most isolated galaxies in the Local Group. WLM 
has a heliocentric velocity of -130 km (jJackson et al.l 
l2004f) and a modest velocity with respect to the Local 
Group barycentre, —29 km s~^, implying that it may 
have recently passed apocentre and is turning around. 
Additionally, WLM lies out of the Galactic plane, which 
minimizes foreground contamination and reddening. 

High dispersion spectra have been taken for a 
few bright A an d B-type supergiant stars in WLM 
(I Venn et al.l 120031 : iBresolin et all l2006t lUrbaneia etHI 
|2008| ). Detailed analyses of these stars provide 
the present day met allicities and ab undance ratios 

([Fe/H]= -0.38 ± .2, iVenn et al.l[200l [O/H] 0.85, 

Bresolinet all 120061 : [Z]= -0.87 ± 0.06. lUrbaneja et all 



20081 ) , but offer little information on the intermediate or 



old age populations. The older red giant branch stars in 
these isolated dirrs are too faint for high dispersion spec- 
troscopic analyses, even with 8m class telescopes. Stud- 
ies of HII regions from emission line spectroscopy yield 
[0/H | = -0.83 (iSkillman et all [1981 iHodee fc Miiie3 
ll995tlLee et al.ll2005[ ). but provide no information about 
the chemistry of the gas in the early stages of the 
galaxy. Interestingly, WLM has revealed minor dis- 
crepancies in chemistry between the young stars and 
H II regions - possibly due to inhomogeneous mixing 
(e.g., IVenn et al. I [20031 foun d [O/H] = -0.21 in one A 
supergiant: iLee et al.ll2005l find [O/H] = -0.83 ± 0.06 
from emission line measurements in two HII regions, 
which included the [OIII] A 4363 line). Neutral gas 
studies in WLM have been used to map the HI enve- 
lope e xtent and small scale spatial and velocity struc- 
tures (iHuchtmeier et all 119811: iBarnes fc de BloH |2004 
iJackson et al.ll2004l : lKeplev et al.ll2007D . 

The old p o pulati on in WLM was first sampled by 
iHodge et al.l ( 19991 ) who derived a metallicity from 
isochrone fitting to deep HST imaging of the lone globu- 
lar cluste r WLM-1. HST photome t ry and wide field INT 
imaging jMinniti & Zijlstral [1993; iReikuba et all l2000l : 
iMcComiac hic et al. 200^ also identified young and old 
stellar populations in the form of an extended blue main 
sequence, and a horizontal branch on the GMD. These 
photometric surveys were also used to find the distance 
and reddening to WLM, and to estimate the range in 
metallicity on the RG B. Photometric analy s es of G and 
M stars in WLM by 'BattincU i fc Demersl ([2004) have 
argued against the presence of an ol d extended halo in 
WLM, opposite to the conclusion from lMinniti fc Ziilstral 
(|1997f ) based on the interpretation of their GMDs. How- 
ever, differential reddening within WLM may be affect- 
ing all of these photometric analyse s: the recent Spitzer 
IRAC survey of AGB stars in WLM (jJackson et al.ll2007f) 
has shown the patchy presence of dust throughout WLM. 



The use of the empirically calibrated near infrared cal- 
cium triplet lines provides a new method for studying the 
stellar population in WLM. Situated at A ~8498, 8542, 
8662 A, they are optimally located with minimal contam- 
ination from other spectral features and near the peak in 
flux for these evolved red stars. The summed equiva- 
lent widths (EWs) of these lines are well calibrated to 
allow a representative placement of a star onto a given 
[Fe/H]^ scale, and sensitive enough from medium reso- 
lution spectra to perform well out to large distances in 
the local volume. The metallicity index is also well cor- 
rel ated with the iron abund ances ([Fe/H]) determined 
bv lCarretta fc GrattonI (|1997| ) from high dispersion spec- 
troscopy of Galactic globular clusters. This scale shows 
a linear correl ation with the CaT W' index, unlike the 
Zinn & West (jZinn fc Westlll984l ) scale based upon the 
Q39 spectrophot ometric index. Previous large scale cal- 
ibration studies ("Rutledge et"al]|1997t ICole et al.| [200l 
iGarrera et al.l20 07b: Battag ha et al.ll2007l ) confirmed the 
robustness of the CaT method over a range of ages 
(0.25 < Gyr < 13) and over the metallicity range ex- 
pected for the stars in a gas rich dIrr galaxy (—2.5 < 
[Fe/H] < -1-0.47). A growing number of CaT studies have 
iDcen carried out for several Local Group galaxies, includ- 
ing the Magellanic clouds (Pont et al. 2004; Cole et aTl 
[2005: Grocholski et al. 2006; Carrcra et al. 2007 a) and 
dSph galaxies (iTolstov et al. 2004; B attaglia et al.ll2006l : 
iKoch et al.ll2006l : iBosler et al.,.2007 '). which further tests 
the robustness of the CaT method in different environ- 
ments. 

In the following sections, we discuss the observations, 
the data reduction methods, and spectral analysis of CaT 
spectroscopy adopted in this paper. These sections are 
particularly important since this work represents a CaT 
analysis of some of the faintest RGB stars for which ve- 
locities and metallicities have been determined at mod- 
erate S/N. The main challenge has been to minimize the 
velocity and metallicity errors. Although the final un- 
certainties are slightly larger than CaT surveys of closer 
galaxies, we have been able to determine a new metallic- 
ity distribution function for WLM and characterize the 
spatial and velocity variations in its stellar populations 
for the first time. 

2. OBSERVATIONS AND DATA REDUCTIONS 

2.1. Target Selection and Photometric Calibration 

Relatively isolated and bright stars near the tip of the 
RGB {Vmag ~ 23) wcrc selected from 150 second V 
and I band preimaging expos ures from F0RS2 on the 
VLT (jAppenzeller et al.l [1991 ). The data were bias cor- 
rected and fiatfielded in IRAF'^ and instrumental pho- 
tometry was obtained using DAOPHOT/ALLFRAME 
(|Stetsonlll994l ). Stars were selected from the tip of the 
RGB (TRGB) to those within half a magnitude below 
the RGB tip (instrumental I band), and colors consis- 
tent with RGB membership. This meant that stars were 
selected to have instrumental colours spanning the ap- 
parent width of the RGB (i.e., constrained to a colour 

2 The notation [Fe/H] = log(Fe/H), - log(Fe/H)Q 
^ IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatory, which is operated 
by the Association of Universities for Research in Astronomy, Inc. , 
under cooperative agreement with the National Science Foundation 
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range of 1.1 magnitudes) to avoid heavy contamination 
by red supergiants and M stars. Our targets were se- 
lected to encompass a broad area of the galaxy's high 
density (gas and stellar) regions, as well as lower density 
outer areas. 

The F0RS2 preimaging photometry was then matched 
to extant INT WFC V, i band photometry of a 
0.25 square degre e regio n centered on WLM (see 
iMcConnachie et all (|2005l ) for further details), due to 
the preimaging being taken under non-photometric con- 
ditions. This allowed us to use standard magnitudes and 
colours for our analysis. A list of INT WFC stars for 
matching was created based on colour, magnitude, and 
location, and these stars were further filtered based on 
sharpness and ellipse parameters. As the stars were uni- 
formly selected to be bright and isolated, we acknowl- 
edge that a selection effect may appear in areas of low 
stellar density. However, this is unavoidable in resolved 
stellar spectroscopy at such large distances, and simply 
necessitates that care be taken in the interpretation of 
the results. Four globular clusters (47Tuc, NGC1851, 
NGC1904, M15) spanning a range of metallicities, were 
also chosen for calibration purposes, as described in 
§3.1.1. 

2.2. Data Acquisition and Reduction 

The observations for this study of WLM and the four 
calibrating clusters were taken during several nights at 
the VLT in late 2003 (see Table 1). The MXU (Mask ex- 
change Unit) mode was used with the F0RS2 instrument 
at the Cassegrain unit of VLT's UT4 (Yepun) telescope, 
note this was before the instrument was relocated to UTl 
(Antu) in 2004. The 83 RGB steUar targets ranged in 
magnitude from 22.1 < Ymag < 24.0, requiring exposure 
times on the order of 40 minutes for individual images 
even with an 8m class telescope. Slit acquisition images 
were taken for approximately 150 seconds in Bessell I 
band, to confirm the pointing and slit alignment prior 
to science exposures. The parameters for the object sci- 
ence exposures on the observational dates in question are 
shown in Table 2. 

Figure [1] shows a representative view of the 
F0RS2/MXU fields used in this study of WLM. For ref- 
erence, the INT WFC fields are the four red rectangles. 
The F0RS2 instrument was used with the MXU in order 
to provide selectable custom cut slit plates; a configura- 
tion which allowed the preselected RGB stars to be fit 
within 1 X 8 slits to minimize the required observa- 
tion time spent for the program. A spectroscopic order 
separation filter (OG590-I-32) and the standard resolu- 
tion collimator were used in conjunction with a volume 
phased holographic grism to obtain the stellar spectra for 
each slit target. The grism, designated GRIS_1028-|-29, 
provided wavelength coverage from roug hly 7730A to 
9480A with a central wavelength of Ac — 86OOA, and 
a dispersion of 28.3A/ mm. F0RS2 is equipped with a 
mosaic of two red-optimized 2kx4k MIT CCDs (15/xm 
pixels) with very low fringe amplitude in the spectral 
range of the CaT lines. With 2x2 binning and lOOkhz 
readout characteristics, and the above mentioned opti- 
cal path, the effective field of view across the instrument 
was 6.8 X 5.7 . The two component chips have a pixel 
scale of 0.252 /pixel, where the pixel size is 0.86A (2x2 



binning). This setup produced observations with resolv- 
ing power R ~ 3400. The gain for both chips was 0.7 
ADUs/e~, with the readnoise being 2.7e~ and 3.15e~, 
respectively. Each of the large blue boxes in Figure [T] 
is representative of the full two chip CCD field of view. 
The representative seeing conditions for the observations 
at VLT ranged from O.6I" < FWHM < 1.52" over the 
course of our science exposures. 

Once obtained in service mode at the telescope, the 
data were reduced using a variety of standard IRAF 
tasks for bias, overscan and flatfield corrections to the 
two dimensional images. Typically each science image 
used a bias correction combined from 10 individual bias 
exposures, along with 5 individual screen flat field ex- 
posures for each of the chips, which were reduced in- 
dependently. The objects were cleaned for cosmic rays 
using several iterations of the COSMICRAYS task. In or- 
der to minimize any error in the final aperture extrac- 
tion or wavelength calibration, custom tasks were run to 
remove warping in the stellar trace and sky lines (par- 
ticularly the far corners of the chips) due to light path 
distortions or chip alignment issues. A custom IRAF 
script was used to provide a preliminary correction to 
the drooping artifact of the stellar trace as a function of 
row position on the CCD. A second script was then used 
to linearize and orthogonalize the dispersion coordinate 
and spati al coordinate to optim ize subtraction of the sky 
lines (see lGrocholski et al.ll2006l ). Following these correc- 
tions, aperture extraction was run on all science spectra, 
resulting in one dimensional spectral images for all the 
stars in the sample. We used a standard sky line atlas in 
conjunction with the identify and REIDENTIFY tasks to 
provide dispersion solutions for the wavelength calibra- 
tion of the spectra. Typical line fits involved ^ 40 OH 
sky emission features taken from lOsterbroc k fc Martell 
(|I992). The typical accuracy of the wavelength solution 
is within ~ 0.04A based on the rms of the solutions. Fi- 
nally, the task continuum removed the curvature of the 
spectrum and normalized the flux distribution to allow 
for the most accurate measurements of the regions of in- 
terest (the calcium triplet lines). 

Prior to combination, we adjusted the wavelength ze- 
ropoint of each spectrum for the heliocentric correction. 
This was a non-trivial concern for the radial velocity mea- 
surements as our individual exposures were spread over 
several months (See Table 2). A date dependent shift 
was applied to each spectrum, such that their observed 
velocities were corrected back to a common epoch at the 
date of the first observation. This resulted in at most a 
~ 1.2 A correction to the spectra. The final spectra in- 
clude the median of eight individual exposures for each of 
79 RGB stars (four stars were thrown out due to spectral 
contamination or photometry matching issues). The typ- 
ical signal to noise ratio for our combined images ranged 
between 16 < (S/N) < 23 per pixel (1 pixel = 0.86A). 

3. SPECTRAL ANALYSIS 

3.1. Equivalent Width Measurements 

The choice of techniques for equivalent width measure- 
ments had to be carefully considered, given the moder- 
ate signa l to noise o f our c ombined spectra (~ 20). As 
noted bv ICole et al.l ()2004f ). the line wings are typically 
underestimated with a pure Gaussian profile fit, and are 



4 



Leaman et al. 



much more accurately modeled with the sum of a Gaus- 
sian and Lorentzian fit. However, for low signal-to-noise 
spectra, when the line FWHM are on the order of the 
spectral resolution, the contaminating noise features ef- 
fectively nullify any difference between using a Gaussian, 
Lorentzian, or the sum of the two. Therefore an accurate 
characterization of the spectrographic setup is necessary 
in low resolution Calcium triplet spectroscopic studies. 
For this reason, the choice was made to use a simple pixel 
to pixel integration of the line profiles. Due to the low 
signal-to-noise of the spectra, the integration was taken 
over the range where the line wings intersected the global 
continuum level. The continuum normalized spectra had 
an average continuum at unity already, no adjustment 
was made to the flux zeropoint. The integration yields a 
wavelength, core density, and equivalent width for each 
line. Multiple independent measurements were made by 
remeasuring the same line repeatedly, with a systematic 
deviation in equivalent width of approximately 1%. Fig- 
ure [5] shows a sample spectrum in the region of the Ca 
II triplet lines for one star in our sample. 

As a consistency check, we performed simultaneous 
measurements on the reduced calibrating clusters (which 
had a much higher signal to noise) using both the pixel 
to pixel integration method and the profile fitting pro- 
grams from Cole et al. (2004). The linear fit between 
the methods' measurements (see Figure [3]) was used to 
develop a transformation function wh ich put our equiva- 
lent widths onto the same scale as in iCole et al.l (|200l . 
thereby allowing us to adopt a variety of their calibra- 
tions as described in the next section. 

3.1.1. Placement onto the Metallicity Scale 

lArmandroff fc Da Costal (|1991l ) first showed the use- 
fulness of the CaT feature as a metallicity indicator in 
individual evolved giant branch stars. Using the summed 
equivalent widths of the CaT lines, as well as the star's V 
magnitude above the horizontal branch, a reduced equiv- 
alent width could be formulated. When observed for 
Galactic globular cluster stars, which also had high dis- 
persion spectroscopic [Fe/H] estimates, this low resolu- 
tion CaT estimator could be used to prov ide an empirical 
meta l licity index Sev eral studies (e.g., iRutledge et al.l 
I1997I: ICole et al.l[20nl have since explored these empir- 
ical calibratio ns over a l arg e range of stellar ages and 
metallicities. ICole et all ([2004) found minimal deterio- 
ration in the CaT calibration over extremes in age and 
metallicity. This is particularly relevant for dwarf galaxy 
studies where we expec t a mixed stella r population based 
on CMD analyses fe.g.. lDolphinll2000f ) In this pa, p er, we 
adopt the calibrations determined by ICole et al] (j2004l ) 
based on many clusters which were made with the same 
instrument (F0RS2) at the VLT, and adopted similar 
data reduction techniques as used here. The spread 
in clu ster ages in the full sample studied bv ICole et al.l 
(|2004f ) is ideal for an alyzing a dir r galaxy lik e WLM . 
Rec entlv. iBattagfia et ai.! |2007l) . ICarrera et all (|2007aD . 
and iBattaglia et al.l (|2008f ) have confirmed the robust- 
ness of the CaT spectroscopic method as a means for 
metallicity estimates in dwarf galaxies. These studies 
specifically checked the appropriateness of the CaT index 
as a proxy for [Fe/H] in cases of varying [g/Fe], [Ca/Fel , 
and age. Additiona l ly, stu dies by iRutledge et al.l ()1997D 
and ICenarro et al.1 (|2001[ ) showed that the CaT index 



transforms in a well understood way between different 
authors' studies. 

To use any of the calibrations requires a summed equiv- 
alent width determination for each star. Each of the 
three calcium triplet line measurements were combined 
in an unweighted fashion to yield a summed equivalent 
width per star. 



Y.W = W8498 + 1^8542 + 1^8662 (1) 

The justification to use all three lines will be discussed 
in §3.2.2 and §3.2.3 with respect to our errors. With this 
relation it is now possible to form the calcium index, W' 
defined as: 



W = Y.W + (3{V -Vhb) (2) 

The term in the parentheses provides a correction for 
the changes in Tgjy and log{g) for stars in different 
phases on the red giant branch. A cooler temperature 
and lower surface gravity play non-trivial roles in the 
formation of the CaT line profiles and the continuum 
in these evolutionary stages. Theoretical and ern piri- 
cal work (jJorgensen et al.iri992l : ICenarro et al.ll200l ) has 
confirmed this complicated interplay of the calcium line 
strengths with stellar parameters such as temperature 
and gravity. Thus, this term is important in removing the 
gravity dependence of the lines with respect to the con- 
tinuum in the CaT analysis. Our V magnitudes are taken 
from the INT WFC catalogue and associated database as 
mentioned in §2.1. We ad opted the horizonta l branch at 
Yhb ^ 25.71 ± 0.09 mag (|Reikuba et "aL 2000), and take 
f3 ^ 0.73 ± 0.04A mag- ^ froin ICole et al.. (.2004 ). Using 
the ICarretta fc GrattonI (|1997f) scale, the calcium index 
is converted to a metallicity ([Fe/HjcGg?) as follows: 

[Fe/H]cG97 = (0.362 ± 0.014)1^' - (2.966 ± 0.032) (3) 

The zero poi nt and slope were determined by 
I Cole et aD ()2004l ). and estimated as accurate to < 4%. 
This means that the majority of our uncertainties are 
from other factors, which we will discuss below in the 
following section. 

Figure m shows a plot of the summed calcium II equiv- 
alent widths (SW) against the V magnitude above the 
horizontal branch per star. The fiducial solid lines show 
the spacing in rn etallicity given by the calibration from 
ICole et al.l (I2OOI . 

3.1.2. Error Analysis 

To examine the various sources of uncertainty in our 
calculations, all errors from the equivalent width mea- 
surements to placement on the [Fe/H]cG97 scale were 
recorded. As no standalone error estimate was provided 
by the pixel integration routine used in measuring a line 
equiva lent width, \V , we adopted the Cayrel fo rmula (see 
Eq. 7: ICavre]|[l988l also iBattagha et al.|[2007l) as: 

{AW^)i ~ {1.6{FWHMJx)ie) + (0.10Ty„) (4) 
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where the pixel size is dx = 0.86A, and the average rms 
continuum accuracy is e ~ (s/n) — ■ pixel in- 

tegration measurements did not contain information on 
the FWHM of the line, that value was determined as 
follows. A test set of Gaussian and Lorentzian fits were 
made to the strongest line (8542 A) in 18 stars, using the 
line profile fitting routine SPLOT. A Gaussian-Lorentzian 
blended estimate of the FWHM was estimated from an 
arithmetic average of the two fit values. These were then 
plotted as a function of measured equivalent width (from 
the pixel integration method). A linear regression was 
fit; an estimate of the FWHM to be used in this formula 
for the rest of the stars was then determined based on 
their equivalent width values. The FWHM values for the 
sub-sample of stars used in this step showed a standard 
deviation of afwhm = 1.06A. 

With a full set of equivalent widths and representa- 
tive FWHM values. Equation 4 could be used to find a 
realistic uncertainty in a given line. These steps were 
repeated for all three CaT lines, then added in quadra- 
ture for the error in the summed equivalent width, EW. 
An additional term was added to this in quadrature, to 
account for the variation in equivalent width measuring 
methods used. This term was determined from the rms 
dispersion about the best fit in the transformation to the 
iCole et all (|2004f) equivalent width measurement system 
(see §3.1.1) 

The final expression for the uncertainty in metallicity 
is now formulated into one equation, allowing for simple 
partial derivative based error propagation. Defining an 
equation for metallicity based on the following equation 
of observables or calibrated variables: 

/ = [FelH]cG97 = ci [HW + I3{Vhb ~ V)] ~ (5) 
allowed us to express the total uncertainty as: 

A[Fe/Hr = I (^)' (Aci)^ + (^)' (ASW^)^ 

+ (i)'(A/3)2 + ...} (6) 

The simple mathematical propagation is reasonable and 
justified, as the the variables involved in Equation (6) are 
uncorrelated and independent. This propagation and er- 
ror accounting results in an average uncertainty in metal- 
licity of A[Fe/H] = ±0.25 dex. 

3.1.3. Three Line Justification 

We performed a parallel reduction omitting the weak- 
est line in the calcium series (8498 A). The purpose 
of this exercise was to test whether the uncertainties in 
constructing the calcium index, W , were significantly 
reduced with the rejection of the lowest signal source. 
Composite two-line indices were created in order to com- 
pare the relative errors that occurred for exclusion of 
each line in the creation of the calcium index. The av- 
erage relative error for the full 80 star sample on any of 
the line pairs, ^^""^ where n and m are the first, second 
or third lines of the Ca II triplet, were calculated to be 
= 0.10, 4r^ = 0.12, and ^ 0.11 for the 



pairs. This minimal deviation, suggests that inclusion of 
all three lines is warranted, because the equivalent width 
measures do not dominate the random error budget, and 
thus there is no benefit to dropping the 8498 A line. 

3.2. Radial Velocity Measurements 

Radial velocities were measured from the strong cal- 
cium lines which had previously had a wavelength dis- 
persion solution applied from the sky OH lines. The low 
signal to noise of the individual frames necessitated that 
we perform the cross correlation radial velocity calcula- 
tions on the combined spectra, rather than each individ- 
ual image. As such, heliocentric velocity corrections were 
tailored to the individual exposures and applied prior to 
combining the spectra, due to the long temporal baseline 
(roughly four months) of the observations. Once shifted 
and combined, the spectra were ready for radial velocity 
computation with the aid of template stars and a Fourier 
cross correlation routine (fxcor). A total of 23 tem- 
plate radial velocity stars, observed with the same instru- 
ment setup, were used with the cross correlation routine. 
This computation provided error analysis automatically, 
with the median error in the heliocentric velocities being 
{SVhei) — ±6 km s""'^. Systematic velocity errors due to 
a star's position in the slit were removed by centroid- 
ing the stars relative to the slit centre. The fact that 
this procedure was done on combined spectra resulted 
in small absolute corrections, as the -y/n statistics meant 
that the individual slit errors were minimized in the com- 
bination and correction steps. The typical shift for the 
slit error on an individual exposure is approximately 6-9 
km s^^, and we note that this shift produces negligible 
uncertainties in the equivalent width error. The final av- 
erage absolute corrections to the slit centering errors on 
the combined spectra were on the order of < 1.5 km s~^. 

Our sample of RGB stars has a mean velocity of 
{Vhei) = —130 ± 1 km s~^, identical to the heliocen- 
tric velocity f o r WL M derived from neutr al HI studies 
iJackson et all (|2004[ ): iKeplev et all ()2007f ). Only one 
foreground star was found, with radial velocity Vuei = 
-1-48 km s~^, leaving us with 78 stars total. Given the 
location of WLM with respect to the Galactic plane and 
the colour and magnitude cuts in our preselection rou- 
tine, it is not surprising that there were so few foreground 
objects. 

3.3. Age Derivations 

An estimate of the relative ages of our sample of RGB 
stars was deter mined using the Victoria -Regina stellar 
evolution tracks (VandenBerg et al. 2006). Given the ex- 
tended SFH of WLM (Mateo , 1998: Dolphin ,2000} and 
the gas rich nature of dirr galaxies, we expect a signifi- 
cant range in red giant ages in the dataset. 

A CMD of the target stars in two metallicity bins, 
along with the fiducial sequences for M68 and 47Tuc, 
is shown in Figure [5l As the stars were selected across a 
wide range in (V-I) colours, the spectroscopic CaT [Fe/H] 
estimates allow us to break the age-metallicity degener- 
acy and provide a more complete understanding of the 
evolved populations in WLM. The age derivation proce- 
dure and error estimat es of the relative stellar ages are 
more fully described bv iLeama 3 (|2008D '^. however a brief 

^ This thesis work is freely available for download at: 
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outline is given below. 

To determine an age from the Victoria-Regina mod- 
els, a value of [Fe/H]caT, V, and (V-I) is required for 
each star, as well as global values for distance and red- 
dening to WLM (see Table 1) and an assumption about 
the a-element enhancement level. The [a/Fe] value to 
use for the WLM sample is not clear especially given 
the age and metallicity range of the evolved populations 
here (e.g., dSph stars show [a/Fe]= -1-0.3 for metal-poor 
sta rs, but show subsqlar ra tios for increasing metallic- 
ity, [Siitox!iES[2Mil^nL&_ Hill 2^ [a/Fe]= +0.3 
was adopted based on the results from galactic globular 
clusters which were used to calibrate the CaT method 
here; however the effect of computing the ages with 
[Q;/Fe]= 0.0 was found to be minor. Systematic ef- 
fects (e.g., an evolutionary model's treatment of mixing 
length, convective overshooting, CNO abundances, ra- 
diative diffusion, gravitational settling etc.) between the 
models and observations were removed by comparing the 
isochrones to homogeneous globular cluster photometry 
from the CADC website^. The location of each star in 
the dereddened plane was used to interpolate a best age 
estimate amongst tracks of the appropriate metallicity 
for that star. 

RGB stars can only evolve to luminosities approaching 
the tip of the RGB when the core becomes significantly 
electron-degenerate before the onset of core helium burn- 
ing. Thus there is a lower limit to the age of the first- 
ascent RGB stars that could possibly be included in our 
sample, typically around 1.6 Gyr. At a given metallicity, 
this age marks the blue edge of the colour distribution 
of RGB stars. In cases where the star was bluer than 
allowed by this condition, an arbitrary age marker of 1 
Gyr was assigned as an aid to interpreting the sample 
properties. Similarly, those redward of the oldest avail- 
able track (18 Gyr) were assigned an age of 18 Gyr and 
a significant error (discussed below). As seen in Table 
3, roughly one third of the stars are below the 1.6 Gyr 
cutoff, and most of those are found in the bar field of 
WLM. 

Formal errors on the ages were adopted by folding in 
the uncertainty in the spectroscopic metallicities, as well 
as the uncertainties in the photometry. While other as- 
trophysical factors could impact the positions of the stars 
on the CMD (e.g., variations in relative internal redden- 
ing), we are unable to expli citly quantify these effects. 
Figure 13 of ICole et cd] (|2005) illustrates the effect that 
some of these factors can have on the ages. Given those 
dependencies and the value of a typical metallicity error 
of ±0.25 dex, the random error in age is ^ ±50%. 

The technique used here for deriving ages can be ap- 
plied successfully in globular cluster populations in the 
MW (|Tolstovl [20031) . However studies using identical 
methods in dwarf galaxies have fou nd it much more diffi - 
cult to assign ages to evolved stars. iTolstoy et al.l (|2003l ) 
first showed that the evolved stars in dSph galaxies of- 
ten lay blueward or redward of the full age parameter 
space (see the discussion in §3.2 Tolsto v et al..,2003: also 
lTolstovll2003D. This pheno menon has also been seen in 
the LMC (|Cole et al.ll200"5[ ). This wide range of colours 
is clear in our Figured! where we overlay GGC fiducial 

|http://hdl.handle.net/1828/1325| 



sequences (M68; [Fe/H]~-2.0, and 47 Tuc [Fe/H]— 0.7) 
onto a metallicity binned colour magnitude diagram of 
WLM. Notable are both metal poor and metal rich WLM 
stars^ blueward of the fiducial sequence for M68. The 
majority of these stars are more metal rich than M68, 
suggesting that we are seeing a range of ages in our sam- 
ple younger than typical globular cluster ages. In Figure 
El we also point out one star (STARID 19203), at (V- 
I)o ^ 2.5, and note that a greater uncertainty should be 
applied to this star's parameters as it may be subject to 
excessive reddening, or spectral contamination (notably, 
from weak TiO bands). 

The question that we are left with is whether the 
bluest stars in our age derivations (those with inferred 
ages <1.6 Gyr), are really intermediate-mass stars on 
the early-asymptotic giant branch or are true low-mass 
stars on the first-ascent RGB, influenced by uncertain- 
ties in composition or differential reddening. If the ma- 
jority of the bluest stars are truly 1 Gyr or less, then 
it is reasonable to look for concentrations of their more- 
evolved descendants (e.g., thermally-pulsing AGB stars). 
Among Local Group dwarfs, WLM has been found to 
have an extraordi narily large ratio of carbon - to M-stars, 
C/M= 12.4±3.7 (jBattinelli fc Demersl[200l . This large 
population of carbon stars, along with the high recent 
star formation of WLM, does suggest that the majority of 
stars in our youngest age bin are truly young. Note that 
these young stars are not likely to be differentially red- 
dened AGB or RSG contaminants, as cross-correlation 
with UKIRT WFCAM JHK photometry rules this out. 
However, early-AGB stars ascending from the horizon- 
tal branch towards the thermally-pulsing (carbon star) 
phase can occupy the same CMD space as a bona fide 
RGB star, and this becomes quite likely in a composite 
population such as WLM's. We note that if an RGB star 
is in fact an early AGB star, at a given metallicity the 
age of that star would be ^ 30% older than we infer here 
(for further discussion see ICole et al.ll200l . 

4. ANALYSIS AND DISCUSSION OF WLM 

In this paper, we have determined the [Fc/H] and 
radial velocity values for 78 individual RGB stars in 
WLM. These can be used to examine the structure, kine- 
matics, and chemical evolution history of this galaxy 
for the first time in conjunction with derived age es- 
timates. This i s unique, as all p r evious studies of su- 
pergiant stars (IVenn et all l2004l : iBresolin et all l2006l: 
Urbaneia et all l2008l) and HII regions (|Skillman et al.l 



'l98aiHodge fc Millej[l995tlLee et al.ll2005D sampled the 
young population and offer little insight into the ear- 
lier epochs of formation and evolution of this galaxy . 
Simil arly, ph otometric studies (Minniti fc Zii lstralll997t 
,Hodge et al.l [199 9; McConnachie et al, 2005) were only 
able to provide global views of the metal poor population, 
and were subject to degeneracies in age and metallicity. 

To begin our analysis, elliptical radii were determined 
for all of the stars in our sample. The position angle 
of WLM was taken t o be 181° (Jackson et al . 2004) and 
the eccentricity 0.59 ([Abies fc AbledligTTlT ^fter trans- 
forming the equatorial coordinates into the transverse 

^ In this, and other representations throughout the paper, the 
metallicity split to characterize "metal rich" and "metal poor" 



url: |http^77www3Tcan^ccn^Tiia-iha.nrc-cnrc.gc.ca/communit y/STE'f5gfflf^Mf^ 



*le at [Fe/H] = -1.27 dex 
; value of our sample. 
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^, r/ plane, and calculating the axial components of the 
ellipse for a given star, the elliptical radius (reii) was 
determined from the geometric mean of the semi-major 
and semi-minor axes. These elliptical radii have uncer- 
tainties due to the inclination with which WLM appears 
to ob servers (69°, which is not assumed in the calcula- 
tions; I Abies fc AblesllTOTTl ). A given star may of course 
sit off the major plane of the galaxy - thus the calculated 
elliptical radii has an intrinsic degeneracy with the star's 
true distance from the plane of the disk. 

A summary of the spatial and photometric properties 
for the 78 RGB stars in this sample are listed in Table 
3. This includes the metallicity, age, radial velocity, and 
elliptical radii values as well. 

4.1. Chemistry 

The metallicity distribution function from our full 
sample of stars provides the first detailed insight into 
the chemistry of the evolved population of WLM. The 
top panel of Figure [6] shows the MDF for the WLM 
sample subdivided into the bar and north fields. The 
global mean (-1.28 ± 0.03, la = 0.37)^ and peak (~ 
— 1.30) values of the MDF are slightly higher than the 
photometric deter minations, but con sistent within er- 
rors. For example. iHodge et al.l (|1999( ) derived [Fe/H] = 
-1.51±0.09 by isochrone f itting to the lone GC in WLM, 
iMcConnachie et all pOOl found [Fe/H] =-1.5 ±0.2 (as- 
suming no alpha-element enhancem ents) from empir- 
ically calibrated RGB colours, and iMinniti fc Ziilstr"al 
(|1997D found [Fe/H] = -1.45±0.2, again from the photo- 
metric properties of the RGB. These photometric metal- 
licity estimates essentially give lower limits, as younger 
stars will be shifted bluer on the RGB. Our higher results 
confirm that effect - as shown by our higher metallicities 
and range of ages for this composite population. 

The MDF in Figure [5] shows a sharp drop off at the low 
end, with very few metal poor stars below [Fe/H] < —1.8. 
iDolphinI (20001) suggested that there was a significant 
star formation event 9-12 Gyr ago with an enrichment 
up to [Fe/H]o = —2.18 ± 0.28, however those stars do 
not appear in our sample. The low metallicity dropoff is 
primarily due to observational difficulties in attempting 
to sample a sizeable relative number of low [Fe/H] stars. 
Metal poor stars at this evolutionary locus are difficult 
to find in any galaxy due to a population bias resulting 
from mass depen dent stellar evolutionary timescales (c.f. 
ICole et al]|2008[ ). For a given time interval, more young, 
high mass, high metallicity stars will evolve to the TRGB 
than older, low mass, metal poor stars in the same time 
interval. The dominant fraction of young stars at the 
TRGB is large enough to outnumber the metal poor old 
stars, despite the contradictory initial relative numbers 
of the IMF. This effect will be enhanced in the presence 
of any age-metallicity relationship and is still present for 
populations with extended SFHs like WLM. In order to 
quantify the probability that the metallicities in our sam- 
ple are drawn from a Gaussian parent distribution (which 
we don't expect due to the extended SFH, but wish 
to quantify), we performed a Kolmogorov-Smirnov (K- 
S) test on the cumulative metallicity distribution func- 

''' It should be noted our observed spread in metallicities should 
be compared to the intrinsic uncertainty in obtaining our individual 
[Fe/H] estimates (i.e., A[Fe/H]=±0.25 dex) 



tion (CMDF). The probability that the distribution was 
Gaussian is ~ 8 x 10^^^. More interesting is a compar- 
ison to a metallicity distribution predicted by a simple 
closed or leaky box chemical evolution model. The bot- 
tom panel of Figure [6] shows our cumulative metallicity 
function compared to predicted distributions from leaky 
box models with a range of values for the parameter 'c', 
which controls the effective yield. The K-S probabilities 
that our distribution is drawn from one of the simple one 
zone models is < 4% for all models shown. Our distri- 
bution clearly shows a lower number of metal poor stars 
compared to these models. 

Not surprisingly, the young stars that make up the 
metal rich half of our sample arc found preferentially in 
the bar field of WLM (see Figure [6] top panel) . Sev- 
eral tracers of high active SF are observed in this in- 
ner field, notably the strong and numerous HII re gions 
(|Hodge fc Milleiifiggstr^ngblood fc Huntedll999D con- 
sistent again with the most enriched material being found 

at low Tell- 

This is more clearly visualized in the top panel of Fig- 
ure [3 where the mean metallicity decreases by ~0.3 dex 
with distance from the galactic core. From that fig- 
ure we find a radial metallicity gradient of ^^"^g/'^^ = 
— 0.14± 0.02 dex/kpc. This smooth gradient may how- 
ever be a line of sight superposition of two (or more) 
spatially localized subpopulations, mimicking a smooth 
radial profile. 

This population change is c onsistent with the result s 
from the photometric study bv lMinniti fc Ziilstral (|1997f ). 
The [Fe/H] values show a dispersion in metallicity that is 
roughly constant over the sampled region of the galaxy 
(cr = 0.32), similar to previous analyses. A two sided 
Kolmogorov-Smirnov test of the bar and north field cu- 
mulative metallicity distributions reveals that they have 
a very small probability of being drawn from the same 
parent distribution {P — 0.005). Splitting our sample 
into those stars with reii < 0.1 deg. and reii > 0.1 deg. 
and running the same K-S test, results in a similar prob- 
ability (P = 0.024). This would appear to strengthen the 
interpretation of distinct spatial populations in WLM, as 
the inner and outer CMDFs show statistically different 
properties - not just different mean metallicities. 

The bottom panel of Figure [7] shows the cumulative 
radial number density (CDF) of metal poor and metal 
rich stars as a function of elliptical radius from the cen- 
ter of WLM. The metal poor stars are found in more 
uniform numbers independent of radial position within 
the galaxy, whereas the more metal enhanced stars are 
found preferentially closer to the center of WLM. 

In Figure [H we compare the unbinned CMDFs for 
three age bins of our sample. In each bin, we have 
compared the WLM subsample to predicted CMDFs 
from a leaky, and closed box chemical evolution model. 
While our sample deviates from the simple models in 
the youngest age bins, the oldest age bin (i > 6 Gyr) 
shows a large probability of being drawn from a distri- 
bution approximated by the closed box model. We might 
naively expect that our oldest stars would have a charac- 
teristic shape to their metallicity distribution that would 
be more closely approximated by a closed box scenario, 
where little infall or outflow occurs. The deviations at 
younger ages presumably occur with the increase of ac- 
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creted gas, and as such, the CMDF of our younger stars 
begins to vary from the shape predicted by a closed box 
model. However, there is an additional astrophysical ef- 
fect that must be considered. As mentioned in para- 
graph two of this section, the evolutionary flux through 
the TRGB results in higher numbers of young stars in 
that locus than older, low mass, metal poor stars. This 
effect means that our 78 star sample is a slightly skewed 
example of the true [Fe/H] distribution at the TRGB. 
As such, one should keep in mind that the probabili- 
ties found in the K-S tests are really describing how our 
sample is comparing to the closed box models - presum- 
ably the true metallicity distribution of WLM would have 
more metal poor stars, and the K-S probabilities would 
be higher. This population bias is much more prominent 
in the youngest age bin, and at the older bin is less of a 
factor, and so these stars may result in a more represen- 
tative sample of WLM's metallicity distribution. With 
larger numbers of stars it will be possible to quantify 
how much of the probability change with age is due to 
an evolutionary population bias, and how much may be 
attributed to infall of pristine gas. 

4.2. Chemodynamics 

Figure [H] shows the results of our radial velocity anal- 
ysis, presented as two velocity distribution functions 
(VDFs) for the bar and north field. Each field has been 
further subdivided into metal poor and metal rich sub- 
samples. We note the peaked, Gaussian distributions 
around the mean values of (vnorth) = —136 ± 3 km 
s~^ and (vbar) = —126 ± 3 km s~^ for the north and 
bar field samples respectively. Comparing these subsam- 
ples shows that the mean velocities of the north and 
bar, metal poor and metal rich, populations are con- 
sistent (i.e., (vNMp) — (vnmb) = —136 ± 3 km s^-'-; 
(vbmp) = (vbmr) = -126 ± 4 km s^^). This indi- 
cates that the general bulk motion of both subgroups of 
stars in a given field holds in accordance with the group 
mean. However the velocity dispersions of the metallicity 
subpopulations, do not show the same uniform charac- 
teristics. The velocity dispersions of the metal poor sub- 
samples for the north and bar fields are larger than the 
metal rich samples in both fields. We find: cri,(jvAfP) — 19 
km s-\ (J^t^NMR) = 12 km s-\ a^t^BMP) = 20 km 
s~^, an d (Tv(BMR) = 14 km s~^. The Sculptor dSph 
galaxy (jTolstoy et al.l[200^ shows similar chemodynamic 
trends. 

Figure [To] again illustrates that the velocity and metal- 
licity of the stars in the WLM sample are changing with 
respect to position in the galaxy. The north field stars 
are more metal poor and have a higher approach veloc- 
ity relative to the bar field stars. This chemodynamic 
signature is consistent with the interpretation of previ- 
ous photometric studies that found a population gradient 
in WLM. However, our results separate the metallicity- 
velocity-position degeneracies, and we have also deter- 
mined ages. 

4.3. Properties of Each Stellar Population and 
Comparisons to the Gas Dynamics 

The sample is divided into three age bins to examine 
the means and dispersions in the metallicities, velocities, 
ages, and elliptical radii of each stellar population, as 



listed in Table 4. Generally, the older stars are more 
metal poor and kinematically hotter than the younger 
stars. It should be stressed however, that the primary 
variations we see in our data are radial velocity and 
metallicity changes. This is confirmed unambiguously 
by a 4-d principal cor nponent analys is, and a k-means 
clustering analysis fsee lLeamanll2008l §4.5). 

Of equal interest is the examination of the properties 
of the stars which do not show the same kinematic traits 
as the gas (from HI studies). The simplest expectation 
is that the metal rich stars would have formed from ma- 
terial that is more recently dynamically coupled to the 
gas, and would have had less time to become kinemati- 
cally hot. This is borne out by the smaller velocity dis- 
persion among the metal rich sub-sample of our dataset 
(see §4.2). 

In addition to this, we can compare the stellar ve- 
locity gradient to that of the HI velocity gradient. 
iJackson et al.l ()2004i ) found an HI velocity gradient con- 
sistent with a rotation curve, and'Kep lev et al.l (|2007f ) de- 
rived a rotation velocity of ~ 30 km s~^ for the gas, with 
a more steeply varying velocity gradient in the north- 
ern approaching side of the galaxy. A direct comparison 
of the ste llar and gas k i nema tics is shown in Figure 1111 
where the lKeplev et al.l ()2007D HI position- velocity curve 
is overlaid on our stellar velocities. The average veloc- 
ity separation between the fiat central part of our stellar 
data and the approaching northern section is ~ 15 km 
s~^. This is statistically significant given our errors of 
±6 km s~^ per star and error on the mean of ± 3 km 
s^^. However the gradient is shallower and only half that 
of the gas rotation velocity. An apparent spatial lag is 
also seen, although the rotation center is consistent with 
the HI studies. These types of stellar and gas velocity 
discrepancies may be common in dwarf irregular galax- 
ies, whether isolated or interacting, as witnessed by the 
massive star kinematics in the Small Magellanic Gloud 
reported by I Evans" fc HowarthI (p008ll 

Clearly the stellar members in this sample are not 
in the same velocity space as the bulk of the rotating 
gaseous component of WLM. In the northern part of the 
galaxy, the stars lag the gas rotation by as much as ~ 35 
km s~^. At the southern most part of our fields, the 
stars show a velocity offset of roughly the same magni- 
tude, but opposite sign. To fully explain the complex 
velocity decoupling seen between the stellar and gaseous 
components of WLM will take a larger sample, but Fig- 
ure [11] nicely illustrates the lack of coherent motion. 

The parameter ^V(^star-gas} has been calculated to 
quantify the difference in velocity of the stars from the 
gaseous components. For each star an offset in velocity 
was found by subtracting the stars' velocity from that of 
the gas velocity at the same right ascension and declina- 
tion in WLM, according to an HI velocity map provided 
by Dr. A. Kepley (priv. comm.). This was done to 
search for subsamples of stars that track the gas velocity 
the closest. If the HI gas in WLM forms a thin rotat- 
ing plane, then the youngest stars in our sample may 
track this better, than say, the oldest stars In Figure 
[T2l ^V(^star-gas) is plottcd against Vhel, age, and metal- 
licity. While there is some substructure in these repre- 
sentations, no clear trends emerge. This suggests that 
the stellar occupation of the HI velocity space is essen- 
tially random - that all stellar age bins / populations have 
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been dynamically excited from the HI bulk motion. This 
would be consistent with the stellar population of WLM 
laying in a thick disk and/or spheroidal configuration. 
Even the youngest stars show a range in ^V(star-gas)j 
thus either the gas and stars were never coupled (e.g., 
if the gas was recently accreted) or the mechanism re- 
sponsible for dynamically decoupling the stellar veloci- 
ties from the gas must have been present within the last 
few Gyr. 

A bar in WLM's recent history could provide a nat- 
ural explanation to the lack of coherent chemistry and 
ages amongst the stars that track the gas velocities the 
closest (recall Figure [TT|) . Bars will dynamically stir a 
stellar population in a low mass galaxy such as WLM, 
and this also offers an explanation for the large velocity 
dispersions and vel ocity offsets in our samp le when com- 
pared to the gas. lOlsen fc MassevI ()2007[ ) noted tidal 
heating may be responsible for a dynamically hot com- 
ponent in the stellar disk of the LMC. However, given 
the isolation of WLM, a past bar may be a more fea- 
sible explanation for the dynamically hot component of 
WLM's stellar body. 

Small scale HI features, includin g a b ar, have been 
invest igated by iKeplev et al.l ()2007t ) and iJackson et al.l 
(|2004f ). Bot h identify a c entral ly evacuated region 
in the disk. iKeplev et al.l ()2007[ ) identify this as a 
ring or "hook" feature related to concentrated star for- 
mation and subsequ ent supernovae explosions, while 
IJackson et al.l ()2004D sug gest either a c entral bar or a 
massive blowout region. IKeplev et al.l (|2007[ ) place an 
upper limit on the age of the ring at 128 Myr, which is 
much younger than the metal rich, young stars in our 
sample (< 1.6 Gyr). However, these young stars also 
show a statistically significant peak in their radial veloc- 
ity dispersion profile at the extent of the ring feature (see 
iLeaman 2008). If these metal rich (and youngest) stars 
are responding dynamically to this feature, then this sug- 
gests the underlying physical mechanism is a resonance 
with longer term effects t han sampled by th e gas fea- 
tures alone. Interestingly, ' Keplev et al.l ()2007D estimate 
that as much as 20% of the total mass of WLM may 
be ass ociated with this gas feature. While IKeplev et al.l 
()2007( ) note the HI isovelocity contours do not show the 
characteristic 's' shape of a bar, nor is the WLM IRAC 
data suggestive of one - the resonances affecting older 
stellar populations set up up by a past bar may outlive 
either of those indicators. 

The ratio of rotation velocity to velocity dispersion 
^itot^ can provide information on the amount that a 
galaxy's structure is rotationally supported. In Table 5, 
we compare this ratio for our full sample of RGB stars, 
and various subpopulations. The stellar velocity ratio 
for the full sa mple is l ower t han that derived for the gas 
kinematics in iMatecl (|1998f l. although the gas velocity 
dispersion and the central density were assumed quanti- 
ties. Our stellar rotational velocities were estimated by 
linearly fitting the heliocentric velocities versus elliptical 
radii. As expected, the metal poor, old, and "halo-like" 
subdivisions produce lower relative ratios than the 
other splits. These results suggest that the older stellar 
population in WLM is more pressure supported, while 
the younger stellar population (and gas) is more rota- 
tionally supported. While all of our stars lay off the 



plane occupied by the gas in position- velocity space, and 
appear dynamically excited and decoupled from the gas 
motion - we see a range of for the stellar subpopu- 
lations within that thick disk of stars itself. 

While all of these results are intriguing, we stress that 
this work is limited by a very small statistical sample, 
particularly our last result since the older stars that sug- 
gest that stellar population is pressure supported rests on 
only 13 individual RGB stars with t > 6 Gyr. A larger 
sample of stars is necessary to examine this and other 
results, and ultimately test the conditions of the isolated 
dwarf galaxies from the earliest epochs of star formation 
to the present. We have shown with this pilot survey that 
spectroscopic analysis of evolved WLM stars is possible, 
and with future observations on 8m class telesecopes we 
will be able to increase the sample size siginifcantly with 
additional fields. 

5. SUMMARY 

We have presented the first detailed analysis of the Ca 
II triplet spectra of 78 individual RGB stars in the iso- 
lated dwarf irregular galaxy WLM. Medium resolution 
spectra taken with the F0RS2 spectrograph at the VLT 
in two fields in this galaxy (a central bar region, and a 
northern field) were coadded to improve the SNR > 20. 
This is the first analysis of such faint RGB stars for both 
velocity and metallicity measurements, and therefore we 
have been especially careful in our data reduction tech- 
niques and error propagation calculations. The equiv- 
alent widths were converted int o metallici t ies for each 
star using the calibration from iGole et all (l2004t ) with 
final errors ~ ±0.25 dex per star. The heliocentric ve- 
locities were determined to within ±6 km (random 
error). Relative ages were derived from theoretical stel- 
lar evolution models, however differential extinction and 
evolutionary effects may bias our age distribution. 

The north field and bar field metallicity distribu- 
tions have mean metallicities and dispersions of [Fe/H] = 
-1.45 ± 0.04, cr = 0.33 from 34 stars, and [Fe/H] = 
-1.14 ± 0.04, a = 0.39 from 44 stars, respectively. Both 
metallicity distributions are narrowly constrained with 
a range of approximately A[Fe/H]~ 1.0. The stars in 
our sample show radial metallicity and velocity gradi- 
ents, but interestingly the stellar rotation velocity is only 
approximately half that of the gas. The stellar velocities 
appear to be kinematically decoupled and in a thicker 
disk than the gas. Detailed comparisons with the HI 
gas, along with the isolation of WLM, suggest that the 
gas and stars were never coupled (e.g., recent HI infall) 
or that the stellar dynamics have been affected by reso- 
nances associated with a past bar. 

WLM shows signatures of a metal poor population 
that is spatially and temporally distinct, and kinemat- 
ically hotter than the young population at all radii - per- 
haps indicative of an extended halo. The velocities of 
the older stars suggest that this evolved stellar popula- 
tion in WLM is pressure supported, as opposed to the 
younger stars and gas disk which appear more rotation- 
ally supported. This latter result rests on a very small 
number of old RGB stars (13 with t > 6 Gyr), but could 
have significant consequences for constraining theories of 
the structure and dynamics of dwarf galaxies in the early 
Universe. A larger sample and survey size is necessary to 
confirm these results and further characterize the stellar 
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populations in WLM. 
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Fig. 1. — Digitized Sky Survey SERC-J image of WLM. The total image is approximately 45' X 45', with North being up and East 
to the left. The relative loeations of the north and bar fields (blue) from this FORS2 spectroscopic work are show n. The fou r red boxes 
indicate the fields imaged by the INT WFC survey (see [McConnachie et al. 2005) Th e two A-type supergian ts from I Venn et al . (2003) are 
located approximate ly in the centre of our bar field, along with the HII regions from lHodge He Miller! 119951 ). and the B supergiants from 
IBresolin et all ||2006|) . 




Fig. 2. — Sample combined spectra for one star (STARID 28531) in the WLM FORS2 north field {top left). The prominent calcium 
absorption features are visible at AA ~ 8498A, 8542A, and 8662 A even in this (S/N)~22 image. The y-axis are relative flux units. A 
dotted line has been drawn at the global continuum level of unity. 
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Fig. 3. — Co mparison of EEW measurements using the pixel integration method described in this work, and the Gaussian + Lorentzian 
profile fits from lCole et alj II2004I ) . for two extreme metal poor and rich calibrating galactic globular clusters. The results of a global linear 
least squares fit to the two cluster stars is shown as the dot dashed line. The rms dispersion about the fit is 0.22 A. 
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Fig. 4. — Summed equivalent width of the Calcium II triple t lines (SW ) versus the V magnitude above the horizontal branch. The solid 
lines show constant metallicity according to the calibration bv lCole et al.l ||2004 ) which we used to derive our final [Fe/H] values (see text). 



Metallicities and Kinematics in WLM 



15 



6 



4 - 



> 



-2 - 



- 



2 



L I I 



"V" 



• ^ ■ " • ..... . . " . . 

.... I' ./.^I':kM.2. v< "'.8"' " 

if » . , 



. »• . . 

». if . " 

I " S' . ' •..,"4 

_ . J .... J 



t a 




■ [Fe/H] < -1.27 



-1 







2 



3 



(V-l) 
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([Fe/H] ~ -0.7) and M68 ([Fe/H]~ -2.0 green dashed line). The full catalogue of INT WFC photometry is shown as the black dots. It can 
be seen that there are both metal poor and metal rich stars blueward of the fiducial sequence for M68. This indicates that wc arc most 
likely sampling a younger population in some areas. More information on age trends can be found in §4.1. 
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Fig. 6. — The metallicity distribution function for the two FORS2/MXU fields in WLM {to-p panel). The bin size is 0.2 dex, and fitted 
Gaussians with parameters shown are indicated by the dashed fines. Note that the mean values for the two fields d ifi'er by ~ . 3 dex , 
and there is an absence of extreme or even moderately metal poor stars in either field. The [Fe/H] estimates from the lVenn et al.l l |2003 ) 
supergiants are shown by the solid arrow. The cumulative metallicity distribution function (unbinned) for our WLM sample is shown in 
the bottom panel. The coloured lines represent cumulative distributions for simple leaky box chemical evolution models, with different 
values (c) controlling the effective yield). Corresponding Kolmogorov-Smirnov D-statistics (D), and probabilities (P) are shown as well. 
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Fig. 7. — [Fe/H] vs. elliptical radii for the full sample of RGB stars in WLM {top panel). A radial change is clearly evident, and 
quantitatively confirms the interpretation of the disk of WLM being a region of higher SF in more recent times relative to the outer galactic 
regions. Cumulative distribution functions vs. elliptical radii for the metal rich and metal poor subsamples {bottom panel). The metal 
poor half of the sample rises more linearly over the radial extent of our survey, compared to the more centrally concentrated metal rich 
subpopulation. Shown are the K-S D statistic, and probability for a two sided K-S test of the two distributions. 
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Fig. 8. — Cumulative metallicity distributions for our WLM sample, subdivided by ages: t < 1 Gyr {top panel), 1 < t < 6 (middle panel), 
t > 6 Gyr (bottom panel). Comparison CDFs based on simple closed (green), and leaky box (blue) chemical evolution models are shown 
for each age set. The two sided Kolmogorov-Smirnov D statistics are shown, as well as the probability that our metallicity distributions are 
drawn from one of the theoretical distributions. The probability rises steeply with age for the closed box model, with P ^ 14%, 46%, 84% 
in the three bins. 
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Fig. 9. — Velocity distributions for metallicity segregated subsamples in the north field {top panel), and the bar field {bottom panel). 
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Fig. 10. — Plot of the heliocentric velocity and metallicity of the 78 RGB stars in our two WLM fields. The mean hai(green triangles) 
and north {open squares) field values are noticeably offset in both chemical and velocity space. The velocity offset between the two fields is 
shown more clearly by the histogram in the subpanel. Errors are ±0.25 dex, and ±6 km respectively per star, and shown for comparison 
on the simulated legend points at the top. 
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Fig. 11. — {top panel) Co mparison of the heliocentric stellar velocity data for the metal rich (red), metal poor (blue) WLM subpopulations. 
Shown in green is the Kc plev et al.l I I2007I) HI gas position-velocity curve. Neither subpopulation tracks the gas particularly well, and the 
velocity gradient is only half that seen for the HI data (see text), (bottom panel) Similar visualization, but with the WLM stars subdivided 
into three age bins: t < 1 Gyr (circles), 1 < t < 6Gyr (triangles), t > 6 Gyr (squares). Negative offsets along the major axis correspond 
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plotted versus [Fe/H] and age. No clear trend exists between the star-gas velocity offset and the metallicity or age values. 
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TABLE 1 
WLM Properties 



Quantity 



Value 



Reference 



il,b) 
E(B-V) 
Distance 
Eccentricity 
Position Angle 
Heliocentric Velocity (vJ^Ji) 
Rotation Velocity (v^J^) 
Mv 



M. 



dyn 



Mhi 



(75.85, -73.63) 
0.035 (mag.) 
932 ± 33 kpc 
0.59 
181 (deg.) 
-130 km s~i 
30 km s^i 
14.1 (mag.) 
2.16 X 10'^ Mq 
(6.3 ±0.3) X 10^ Mq 
-1.45 ±0.2 dex 



Gallouet et al. (1975) 
McConnachie et al. ( 2005) 
^IcConnachie et al. (2005) 
Abie s & Abies ( 1977) ' 
Jack son et al. (200^ 
Jackson et al. (2004) 
Keplev et al. (2007) 
vandenBergh (1994) 
Keplev e t al. (2002) 
Keplev e t al. (2007) 
Minniti fc Ziilstra (199T\ 
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TABLE 2 
Observational Parameters 



Target Mag. Range Distance (kpc) R.A. Dec. Exp. Time (minutes) 



WLM-N 22.1 < Vmag < 24.0 970'' 00:01:58 -15:21:47 320 

WLM-B 22.2 < Vmag < 23.2 970^^ 00:01:58 -15:28:30 320 

NGC 104 11 < Vmag < 14 4.5 00:24:05 -72:04:51 1 

NGC 1851 12 < Vmag < 16 12.1 05:14:06 -40:02:50 1 

NGC 1904 12 < Vmag < 16 12.9 05:24:10 -24:31:27 1 

NGC 7078 12 < Vmag < 16 10.3 21:29:58 +12:10:01 1 



Note. — The location and distance data for the four calibrating globular clusters is taken from the [Harris l ll996f) catalogue which can 
be found online at http://physwww.mcmaster.ca/~harris/mwgc.dat 
From lGieren et al.i ((20081) 
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TABLE 3 

Selected Parameters For WLM Stellar Sample 



STARID R.A. 


Dec. 


reii V 


I Wi 


W2 


W3 


[Fe/H]car 


A[Fe/H] 


V^ei AVhe! Age 


(degrees) 


(degrees) 


(degrees) (mag.) 


(mag.) (A) 


(A) 


(A) 


(dex) 


(dex) 


(km s-i) (km s"!) (Gyr) 



27626 


0.5293 


-15, 


.5334 


0, 


,0978074 


22, 


.54 


21 


.43 


1 


.076 


2.108 


2 


.781 


-1, 


28986 


0.4753 


-15, 


.5315 


0, 


,0889742 


22, 


.82 


21 


.36 


1 


.209 


3.408 


3 


.208 


-0, 


29674 


0.4773 


-15, 


.5292 


0, 


,0857606 


22, 


.96 


21 


.25 


1 


.250 


3.316 


3, 


.527 


-0, 


28080 


0.4978 


-15, 


.5266 


0, 


,0814584 
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0.21 
0.26 
0.27 
0.22 
0.27 
0.23 
0.25 
0.27 
0.28 
0.21 
0.25 
0.25 
0.22 
0.26 
0.20 
0.22 
0.24 
0.22 
0.25 
0.25 
0.26 
0.25 
0.13 
0.25 
0.24 
0.25 
0.25 
0.27 
0.23 
0.23 
0.22 
0.24 
0.25 
0.31 
0.24 
0.25 
0.25 
0.17 
0.26 
0.24 
0.27 
0.25 
0.28 
0.29 
0.23 
0.21 
0.22 
0.23 
0.24 
0.23 
0.24 
0.22 
0.23 
0.21 
0.21 
0.27 
0.24 
0.21 
0.27 
0.26 
0.22 
0.19 
0.24 
0.28 
0.25 
0.23 
0.16 
0.21 
0.22 
0.21 
0.21 
0.21 
0.22 
0.26 
0.20 
0.23 
0.23 
0.24 



-84 
-112 
-145 
-123 
-128 
-111 
-141 
-119 
-128 
-131 
-142 
-147 
-128 
-134 
-161 
-161 
-129 
-98 
-119 
-107 
-120 
-116 
-109 
-123 
-139 
-129 
-124 
-123 
-114 
-133 
-128 
-132 
-111 
-115 
-130 
-145 
-141 
-109 
-130 
-126 
-107 
-98 
-149 
-156 
-116 
-133 
-144 
-148 
-130 
-103 
-115 
-137 
-139 
-118 
-123 
-130 
-153 
-107 
-124 
-127 
-126 
-189 
-117 
-149 
-146 
-144 
-148 
-111 
-140 
-130 
-148 
-162 
-144 
-132 
-1.30 
-151 
-151 
-147 



< 1.6 
2.4 
5.0 
14.0 
2.2 
2.0 

< 1.6 
3.8 

< 1.6 
2.6 
3.0 

< 1.6 
8.0 

< 1.6 
18.0 
3.2 

< 1.6 

< 1.6 
5.0 
4.2 
3.4 

< 1.6 

< 1.6 
4.6 
2.4 

< 1.6 

< 1.6 

< 1.6 
16.0 
2.8 
10.0 

< 1.6 

< 1.6 

< 1.6 

< 1.6 
4.2 
1.8 
18.0 

< 1.6 
2.9 

< 1.6 
1.8 

< 1.6 

< 1.6 
3.5 
18.0 

< 1.6 
14.0 

< 1.6 

< 1.6 
4.2 
18.0 
2.4 
18.0 
1.8 
2.8 

< 1.6 
3.8 
1.8 
2.6 

< 1.6 
18.0 
2.0 
1.8 

< 1.6 

< 1.6 
3.2 
3.2 

< 1.6 

< 1.6 
8.0 
2.6 
2.6 
1.8 
16.0 
3.2 
2.7 

< 1.6 
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TABLE 4 
Binned Age Metallicity Statistics 



Variable 


t < 1 


1 < t < 6 


t > 6 


{[Fe/H]) (dex) 


-1.21 ±0.05 


-1.23 ±0.04 


-1.55 ±0.07 


(^■[Fe/H]) (dex) 


0.42 


0.34 


0.22 


{vhel) (km 


-129 ± 1 


-130 ± 1 


-136 ± 2 


(<T„) (km s-i) 


17 


16 


22 


(Gyr) 


0.2 


3.0 


15 


(reii) (deg.) 


0.07 


0.09 


0.08 


{ar) (deg.) 


0.05 


0.06 


0.05 


n 


30 


35 


13 



Note. — Average values in three age bins, for various parameters of our WLM stellar sample. 
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TABLE 5 

SUBPOPULATION ^^^^ RATIOS 



Component 


Vrof 




Vrot 


Full Sample 


20.78 


17.53 


1.19 


Metal Rich 


20.61 


14.29 


1.44 


Metal Poor 


19.64 


19.91 


0.99 


t < 1 Gyr 


22.55 


17.31 


1.30 


1 < t < 6 Gyr 


17.52 


16.14 


1.09 


t > 6 Gyr 


9.90 


15.20 


0.65 


Hl-like velocities ("disk") 


40.50 


13.50 


3.00 


non-HI-like velocities ("halo") 


22.86 


18.67 


1.23 



Note. — The fi nal two rows separate the sample based on those stars that fall in the same position- velocity space as the 
iKeplev et all l|2007l ') HI velocity data. 



